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CAUSALITY OF PLASMA PERMITTIVITY WITH LOSS AND RESTORING TERMS 


ABSTRACT 


The causality of the complex permittivity for a plasma model, 
in which both loss and restoring terms are included, is shown. 
Discussions as to what occurs to this causality property when the 
loss and restoring terms vanish are presented. Results are shown 
to agree with earlier work. 


INTRODUCTION 


The expression representing the complex permittivity of a 
plasma material is derived through the equation of motion of 
electrons in the presence of an applied electric field [1-2]. In 
many plasma models, the material is assumed inhomogeneous in the 
direction normal to the plane layer and homogeneous in the 
unbounded lateral directions. The inhomogeneity of the material 
is expressed in the complex permittivity representation through 
the plasma and collision frequencies. An electron density which 
is identified with the plasma frequency and a loss term which is 
identified with collision frequency vary as one moves into the 
material. It is this spatial dependency of the plasma and 
collision frequencies which require numerical evaluation of 
Maxwell's equations in the inhomogeneous plasma material as a 
function of the applied electric field frequency. 

Since the complex permittivity is dispersive, its frequency 
behavior for a fixed plasma profile is also of interest. In 
particular, the complex permittivity should be causal; that is, 
the time response for the complex permittivity must exist only for 
times greater than zero. The expressions given in the literature 
representing complex permittivity [1-3], which are derived with 
only losses (collision frequency) assumed in the model, are shown 
to be noncausal [4] . It was shown in reference 4 that the system 
can be made causal by the addition of a term proportional to the 
impulse function 5 (w) . 

The purpose of this paper is to show that the complex 
permittivity representation, in which both loss and restoring 
terms are included in the plasma model, obeys the principle of 
causality. By a proper limiting procedure as the restoring term 
vanishes, the same impulsive term that appears in the solution 
given in reference 4 is recovered. 
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PROOF OF CAUSALITY 


The differential equation which represents the dynamics of 
electron motion in a plasma or dielectric material in the presence 
of an applied electric is given as [5] 
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where m is the mass of an electron, q is the charge of an 
electron, v is the collision frequency (corresponds to loss) of 
the plasma, b is the coefficient of the restoring force term, x(t) 
is the displacement of the electrons, and £(t) is the applied 
electric field. The derivation of the complex permittivity based 
on equation (1) is given in reference 5. Therefore, without 
proof, the complex permittivity for this model is given as 
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where is the plasma frequency, v is the collision frequency of 
the plasma, P = b/m, is the free space permittivity, and u is 
the frequency of the applied electric field. 

Rewrite equation (2) as 
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The inverse Fourier transform is defined as 
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Substituting equation (3) into equation (4), 
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or 


e (t) = e 5 (t ) - 
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The integral in equation (6) can be readily evaluated by contour 
integration, once the poles of its integrand are determined. 

The poles are determined by setting the denominator of the 
integrand in equation (6) to zero. The discriminant obtained by 
using the quadratic formula permits, in general, three solution 

2 2 

sets. The solution set corresponding to the condition v <4(3 is 
considered invalid because as (3 is ultimately allowed to approach 
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zero, v must be real. For the condition v >4(3 , the poles are 
given as 
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which are shown in figure 1. 

jut 

The complex function f(u) = — ^ is meromorphic in the 

u -(3 -jui> 

upper half plane of the complex-u plane and analytic in the lower 
half plane of the complex-u plane. From Cauchy-Goursat theorem, 
the closed contour integral of f(u) for t < 0 is zero [6]; that 
is. 


<j> f(u) du = 0 (9) 

where the closed contour is indicated by the dashed lines in 
figure 1. Rewriting equation (9) over the individual contours 
shown in the figure: 
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It is easily shown that the contribution from the C r< integral is 
zero [6] . Equation (10) now becomes 
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which says that e(t) in equation (6) is zero for t < 0. 

For t > 0 the contour may be closed in the upper half plane 
as indicated by the solid line shown in figure 1. The closed 
contour integral must again be zero; that is. 


<j> f(to) dto = 0 


( 12 ) 


Expanding the total contour into the individual contours as shown 
in figure 1: 
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The contribution from the contour C r> is zero. The contributions 

from second and third integrals are evaluated by residue theory. 
It is easily shown that 
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With these equations, equation (6) becomes 
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for 4/3 2 < i/ 2 and t ^ 0. For 4 /3 2 < i> 2 and t < 0, e(t) equals zero. 

2 2 

The condition 4/3 = v and t > 0 can be evaluated by 
considering the contour integration path as shown in figure 2. In 
this case there exists a pole of order two at co = ji>/2. Evaluating 
by residue theory, 
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and, since there are no poles in the lower half plane, e(t) is zero 
for t < 0. Hence, equation (6) becomes 
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This result can also be obtained from equation (16) through a 
proper limiting procedure. 


Letting (3 approach zero, that is the term corresponding to 
the restoring force vanishes, equation (16) for 4/3 2 < i> 2 becomes 


e(t) 
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which , once normalized with respect to c , is the same as equation 
(7) in reference 4. Letting 13 approach zero in equation (18), 
which implies v goes to zero because v^~ 40 2 , the expression for 
e(t) becomes 


c(t) = e5(t) + cj 2 c t , t £ 0 (20) 

which, once normalized with respect to e , is the same as equation 
(19) in refere nc e 4, 


It is emphasized at this point that e (t) obeys the principle 
of causality whenever both the loss and restoring terms are 
included (see equations (16) and (18)). Allowing /3 to approach 
zero or both 13 and v to approach zero in the general solutions 
(equation (16) and equation (18)), causal solutions are obtained 
(equation (19) and equation (20), respectively). However, if /3=0 
in equation (1), the direct derivation of c (o>) and its inverse 
Fourier transform lead to a noncausal solution (see equations (4) 
and (5) in reference 4). Also, if /3 = 0 and v = 0 in equation (1), 
the direct derivation of c (u>) and its inverse Fourier transform 
lead to a noncausal solution. Reference 4 shows, however, that 
each of these noncausal solutions can be made causal by the 
addition of an impulsive term. The next section shows, through a 
limiting procedure as (3 approaches zero in equation (2), that the 
proper impulsive term given in equation (6) in reference 4 is 
recovered . 
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PERMITTIVITY A£ RESTORING TERM VANISHES 


In this section equation (2) is examined as 13 approaches 
zero. Its inverse Fourier transform is written explicitly as 
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The evaluation of this integral as (3 approaches zero is performed 
by contour integration. The poles of the integrand are 
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where the condition 4/3 < is the only valid solution because v 

is real and greater than zero. These poles are indicated in 
figure 3. 

Denoting the first integral by Ij and the second integral by 
I^ , equation (21) is rewritten as 
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Equation (24) can be evaluated in a straightforward but lengthy 
process using the contour shown in figure 3. For brevity, only 
the results are indicated: 
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(26) 


At this point it might appear that as (3 approaches zero I 2 would 

vanish. This is incorrect because the range of integration 
includes w = 0, and thus the square bracket in equation (25) 
2 4 

varies as 13 /u> and is not zero. Contour integration, using the 
contour shown in figure 3, yields 
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Substituting equations (26) and (27) into equation (23), 
expression e becomes 
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An interesting result can be obtained if equation (25) is 
rewritten as 
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Comparing this equation with equation (27), the following 
definition can be made: 



(30) 


This is one of the many different definitions which can be made 
for impulse functions. Definitions of others are given in 
reference 7. It is observed that as £ approaches zero, therefore, 
equation (2) must become 
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which, once normalized with respect to e, is identical to 
equation (6) in reference 4. 
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CONCLUDING REMARKS 


The causality of the complex permittivity for a plasma model, 
in which both loss and restoring terms are included, has been 
shown. Whenever the loss term or the restoring term or both were 
allowed to vanish in the time representation for this plasma 
model, the solution still obeyed the principle of causality. 
However, direct solutions of the differential equations without 
both of these terms produced noncausal results . Through an 
appropriate limiting procedure for the restoring term vanishing, 
the impulsive term given in an earlier paper was recovered — the 
term required for causality. 
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Figure 1. -Contours for equation(6) with 4|3 2 < v 2 
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Figure 2. -Contours for equation(6) with 4p 2 = t> 2 
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Figure 3. -Contours for equation(21) 
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